A Dynamical Scaling Law for Jet Tomography 
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Medium modifications of parton fragmentation provide a novel tomographic tool for the study of 
the hot and dense matter created in ultrarelativistic nucleus-nucleus collisions. Their quantitative 
analysis, however, is complicated by the strong dynamical expansion of the collision region. Here, 
we establish for the multiple scattering induced gluon radiation spectrum a scaling law which relates 
medium effects in a collision of arbitrary dynamical expansion to that in an equivalent static scenario. 
Based on this scaling, we calculate for typical kinematical values of the RHIC and LHC heavy ion 
programs medium-modified parton fragmentation functions for heavy ion collisions with realistic 
dynamical expansion. 



For the study of the shortlived state of hot and dense 
matter produced in nucleus-nucleus collisions, novel tools 
become available at collider energies. In particular, the 
abundant production of high-pj_ partons in hard pro- 
cesses dominates the high-p^ tails of hadronic particle 
spectra. These partons propagate through the medium 
while fragmenting into hadrons. Basic properties of the 
hot and dense matter such as the average in-medium 
pathlcngth of hard partons and the transverse colour 
field strength (energy density) they experience are thus 
reflected in the medium-dependence of parton fragmen- 
tation. 

Motivated by this idea (l), several groups ||, [|, ||, |f| 

calculated in recent years gluon radiation spectra due to 
medium-induced multiple scattering contributions in or- 
der to understand the medium-dependence of hadronic 
cross sections. The aim of this letter is to extend previ- 
ous studies || |?j to the full gluon radiation spectrum in 
a strongly expanding medium of small finite size, to cal- 
culate for the first time the medium-dependent fragmen- 
tation functions resulting from many soft interactions of 
the hard parton with the medium, and to explore observ- 
able consequences for collider experiments at RHIC and 
LHC. 

Hadronic cross sections for high p± particle production 
are calculated by convoluting the parton distributions 
of the incoming projectiles with the product da h (z,Q 2 ) 
of a perturbatively calculable partonic cross section a q 
and the fragmentation function D h / q (x,Q 2 ) of the pro- 
duced parton, da h (z,Q 2 ) = (^§jf) dy D h / q (x , Q 2 ) dx . 
Here, x = Eh/E q , y = E q /Q and z = Eh/Q denote 
fractions between the virtuality of the hard process Q, 
and the energies of the produced parton and resulting 
hadron. If the produced parton looses with probability 
P(c) an additional fraction e = ^p- of its energy due to 
medium-induced radiation, then the hadronic cross sec- 
tion is given in terms of the medium-modified fragmen- 
tation function [||, |j| 
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The probability that a hard parton looses AE of its initial 
energy due to the independent emission of an arbitrary 
number of n gluons is determined by the medium-induced 
gluon energy spectrum Q 
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We evaluate the quenching weight (||) via its Mellin trans- 
form (sec 1 10 for details), starting from the medium- 
induced BDMPS-Z gluon radiation spectrum 0, pd|, 
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Medium properties enter dl/dio via the product of the 
medium density n(£) of scattering centers times the rel- 
ative strength of the single elastic scattering cross sec- 
tion oc |a (q±)| 2 measured by the dipole cross section 
<r(r) = 2C A /^| ao (qJ| 2 (l_ e ^- r ). We con- 
sider the effect of arbitrarily many soft momentum trans- 
fers to the parton, when the path integral in ( |3|) can 
be evaluated in the saddle point approximation || |J, 
a(r) w Cr 2 . This is complementary to studies which 
consider one additional medium-induced gluon exchange 
via twist-4 matrix elements or up to N < 3 

medium-induced momentum transfers in an opacity ex- 
pansion |, § of (§). 

For the case of a nuclear medium without dynami- 
cal evolution, n(£) = no, the radiation spectrum (|^) 
depends on the partonic in-medium pathlength L and 
the transport coefficient uqC. The latter measures the 
squared average momentum picked up by the partonic 
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projectile per unit pathlength p] . Phenomenological 
estimates O, 14 range typically between noC = 
(50MeV) 2 /fm for normal nuclear matter and n C ~ 
(500MeV) 2 /fm for a hot quark-gluon plasma. Writing 
all energies in units of the characteristic gluon frequency 
uj c = 2ngCL 2 , the probability P(AE/cu c ) depends on 
only one further, dimensionlcss parameter combination 



R = n CL i 



(5) 



We refer to R as "density parameter" since it character- 
izes the initially produced gluon rapidity density (see eq. 
(|ll]) below). 

Fig. Q shows the quenching weights P(AE/lu c ) calcu- 
lated for a static nuclear medium and varying density 
parameters R. In comparison to earlier studies fllCf] , we 
observe as a novel feature the occurence of a discrete 
contribution po in 



P(AE/u c ) = p 5(AE/u c )+p(AE/u c ). 



(6) 



This is a consequence of considering a medium of realistic 
small size and opacity, where the projectile escapes the 
collision region with finite probability po without further 
interaction. With increasing density parameter R, the 
probability po of no interaction tends to zero, and the 
width of p(AE /lj c ) broadens since a larger energy frac- 
tion AE is lost. In [[To|| , an analytical approximation for 
P(e) is given solely for illustrative purposes. It is based 



on the small-w approximation 4-j oc a 



sion valid for large in-medium pathlengths only. We find 
that it captures well the shape of P for large systems 
where P q <C 1. However it shows an unphysical large-e 
tail with infinite first moment J dee P(e). As seen from 
eq. ( flC)| ) discussed below, the accuracy of our study prof- 
its significantly from a better numerical control of this 
large e region. 

The strong longitudinal and - to a lesser extent - trans- 
verse expansion of a heavy ion collision implies a decreas- 
ing density of scattering centers n(£), 



n(OC = n d C 



(7) 



Here, a — characterizes a static medium and a = 1 
corresponds to a one-dimensional, boost-invariant longi- 
tudinal expansion consistent e.g. with hydrodynamical 
simulations. The maximal value rid is reached around 
the formation time £o of the medium which can be set by 
the inverse of the saturation scale p sa t [ fl6| , resulting in « 
0.2 fm/c at RHIC and « 0.1 fm/c at LHC. Since the dif- 
ference between £ and the production time of the hard 
parton is irrelevant for the evaluation of the radiation 
spectrum (|^), we replace the latter in (||) by £o- 

We evaluate the gluon radiation spectrum for dynam- 
ically expanding collision regions over a wide range of 
expansion parameters a e [0 : 1.5] by combining the 
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FIG. 1: The two contributions to the probability (El) that 
a parton looses AE of its energy in the medium: Continu- 
ous part (upper figure) and the discrete probability po in (^) 
that the hard parton escapes the medium without interaction 



analytic solution || of the path-integral (^) for a time- 
dependent harmonic oscillator with the treatment of fi- 
nite size effects j[l)- We observe a simple scaling law 
which relates the gluon radiation spectrum (^|) of a dy- 
namically expanding collision region to an equivalent 
static scenario. The linearly weighed line integral 
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defines the transport coefficient of the equivalent static 
scenario. Fig. || illustrates the accuracy of this scaling 
law for expansion parameters a = 0, 0.5, 1.0, 1.5, and 
different values of the density parameter R — nCL 3 . 
The accuracy improved with increasing density param- 
eter R. On the level of the quenching weight (||) , it was 
better than 10 % in the physically relevant parameter 
range, except for very thin media (R < 100) where en- 
ergy loss effects arc negligible (< 2%). In a subsequent 
publication |l7| ], we shall extend this scaling law to the 
kj_ -differential radiation spectrum, and we shall docu- 
ment a CPU-inexpensive numerical routine which allows 
to implement the quenching weight (Q) for arbitrary val- 
ues of a, ndC and L in event generator studies of hadronic 
spectra. 
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FIG. 2: The medium-induced gluon radiation spectrum lu jr- 
for different dynamical expansion parameters a=0 (solid line), 
a=0.5 (dotted-dashed), a=l (dashed) and a=1.5 (dotted). 
The scaling law (§) is shown for R=400 and R=4000. 



The linear weight in (|8j) implies that scattering centers 
which are further separated from the partonic produc- 
tion point £o are more effective in generating gluon radi- 
ation. For a dynamical expansion following Bjorken scal- 
ing [a = 1 in eq. (Q)] , all timescales contribute equally to 
nC. This indicates that observables related to medium- 
induced gluon radiation give access to the quark-gluon 
plasma lifetime. 

Since the radiation spectrum M- depends on the ex- 
pansion parameter a, the length L and the transport 
coefficient only through nC, dynamical scaling holds au- 
tomatically for the average parton energy loss 



(AE) = JdEEP(E) = Jdcju;^. 



(9) 



This is seen already from the expressions for (AE) de- 
rived for dynamically expanding scenarios by Baier et 
al. [|[ in the dipole approximation and by Gyulassy et 
al. 1 7) for N = 1 scattering center. Moreover, Gyu- 
lassy et al. JtJ conjectured that higher order opacity 
terms (N — 2,3) of (AE) show the same dependence 
on nC . The results presented here go beyond establish- 
ing this conjecture and confirming the result of Ref. ||. 
Our novel finding is that dynamical scaling holds beyond 
the w-integrated average energy loss @ to good numer- 
ical accuracy for the cj-differential gluon radiation spec- 
trum. This is important since a reliable calculation of 
the medium-dependence of hadronic spectra requires ]Tc| ] 
knowledge of the quenching weight (^) beyond the first 
moment (^) [i.e., requires knowledge about More- 
over, the scaling law established here is relevant for the 
angular gluon radiation pattern [ ]l~7[ , since the radiation 
spectrum emitted outside the opening angle <d can be cal- 
culated in the present approach essentially by replacing 
R = n CL 3 -> sin 2 6 n CL 3 . 



We calculate the medium dependence of parton frag- 
mentation functions Df l / q (x,Q 2 ) from (fy) using the 
quenching weights (||) and the LO BKK [fL5| parametriza- 
tion of D h / q (x,Q 2 ). Typical results for the pion frag- 
mentation function of up quarks are shown in Fig. |^. 
For the present study, we identify the virtuality Q of 
Dh/ q ( x ,Q 2 ) w hh the (transverse) initial energy E q of 
the parton. This is justified since E q and Q are of the 
same order, and Df l / q (x,Q 2 ) has a weak logarithmic in- 
dependence while medium-induced effects change as a 
function of e = 4f ~ 0(±). 

As seen in Fig. |3J, the fragmentation function de- 
creases for increasing values of the transport coefficient 
nC, since the probability of a parton of initial energy 
E q to fragment into a hadron of large energy xE q de- 
creases with increasing parton energy loss. The relative 
size of this medium modification changes roughly like 
e = AE ~ O(^j). We note that the calculation based 
on ([j]) is not reliable for small fractions x (x < 0.1 say), 
since it implies for D^J q d \x, Q 2 ) a normalization 



dx x D\"! ed ^ (x) 

h/q \ I 



f dxxD h/q (x) J de(l-e)P(e) 



(10) 

which is a factor J deeP(e) too small. The origin of this 
error is that the hadronized remnants of the medium- 
induced soft radiation are neglected in the definition of 
(|l|). Due to the softness of these remnants, however, the 
true medium modified fragmentation function is under- 
estimated by equation (|l|) for small x only. We expect 
eq. ([!]) to be valid for the calculation of D < ^J q id \x,Q 2 ) 
for x > 0.1. 

To illustrate the effect of medium-modified fragmen- 
tation functions on hadronic cross sections, we plot in 
Fig. ^ the fragmentation function multiplied by x 6 . This 
exploits that hadronic cross sections weigh the fragmen- 

by the initial hard par- 



tation function D^^ d \x,Q 2 



h/q 



tonic cross sections da q /dp 2 L ~ l/p"^^ and thus ef- 
fectively test x n{ ^D ( ™ cd \x,Q 2 ). We choose n = 6 
for illustrative purposes. For interpretation, the position 
of the maximum x max of x 6 D < j^J^ d \x,Q 2 ) corresponds 
to the most likely energy fraction x msx E q of the lead- 
ing hadron. The medium-induced relative reduction of 
x 6 D^^(x, Q 2 ) around its maximum translates into a 
corresponding relative suppression of this contribution 
to the high-px hadronic spectrum at p± ~ x max _E 9 . For 
Q = 10 GeV, e.g., the leading hadron has most likely 
px ~ 5 - 8 GeV, and a reduction of this contribution 
to the pion spectrum by a factor 2 corresponds e.g. to 
L = 7 fm and nC =1-5 fm -3 , see Fig. 0. The transport 
coefficient can be related to the initial gluon rapidity den- 
sity |^, Rj. Uncertainties in this procedure remain to be 
discussed. For the purpose of comparing our parameter 
values directly to results in the recent literature jffl, we 
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FIG. 3: The LO BKK jlj fragmentation function u -> n for 
no medium and the medium-modified fragmentation functions 
for different gluon rapidity densities (see eq. ([H])) and L — 7 
fm. 



use for a Bjorken scaling scenario 



R = nCL 



L 2 dN3 



R\ 



dy 



for a = 1 , 



(11) 



where Ra denotes the nuclear radius. The curves in 
Fig. H thus correspond to nC ks 1, 5 and 10 fm -3 , re- 
spectively. 

A quantitative study of hadronic cross sections and a 
careful discussion of theoretical uncertainties is clearly 
needed for a detailed comparison of the above results to 
measured hadronic spectra. It is nevertheless interesting 
to contrast Fig. Hwith the transverse 7r°-spectrum mea- 
sured at RHIC |l§|]. Taking the above-mentioned val- 
ues at face value, the soft multiple scattering approach 
to parton fragmentation functions relates a medium- 
ind uced suppression factor w 2 at v ~ 5 GeV to an ini- 



tial eluon rapidity density which seems comparable with 



the value 500-1000 extracted in Ref. 0. Thus, while the 
multiple soft BDMS scattering approach without finite 
size treatment was criticised Q| for overestimating the 
effects of parton energy loss below p± w 10 GeV by up 
to an order of magnitude, the present findings and those 
obtained by a lowest order opacity expansion seem com- 
parable. Differences pointed out Q earlier do not allow 
to prefer the N = 1 opacity expansion over the multiple 
soft scattering approach but just indicate that finite size 
effects have to be treated in both approaches with the 
same scrutiny. 

Equation (|ll"| ) suggests that the average transport co- 
efficient nC changes from RHIC to LHC proportional 
to the observed particle multiplicity. Assuming a factor 



~ 5 increase consistent with recent studies JT6|, Fig. || 
suggests the existence of a sweet spot at LHC energies 
in the p^-range around 50 GeV, where theoretical un- 
certainties are much better controlled than at p± < 10 
GeV while suppression factors are still sufficiently large 
(of order 2) to be experimentally accessible. Attaching 
more precise numbers to this discussion lies beyond the 
exploratory calculation presented here. It does not only 
require a quantitative understanding of hadronic spec- 
tra in ultrarelativistic nucleus-nucleus collisions based on 
the knowledge of the underlying partonic cross sections 
and the nuclear parton distribution functions. It also 
requires medium-modified fragmentation functions for a 
realistic distribution of in-medium pathlengths L and 
initial densities. This in turn requires modeling of the 
spatio-temporal evolution of the hot and dense medium. 
The scaling law established here simplifies the inclusion 
of these dynamical effects in quantitative studies. It can 
be extended to the study of the medium-dependence of 
angular radiation patterns jl7| . 

We thank F. Arleo, K. Eskola, A. Kovner, K. Redlich 
and V. Ruuskanen for helpful discussions. C.A.S. is sup- 
ported by a Marie Curie Fellowship no. HPMF-CT-2000- 
01025 of the European Community TMR program. 



1T0 

[11 
[12 

[is; 

[14 



[is; 

[16 

li- 
tis 



M. Gyulassy and X. N. Wang, Nucl. Phys. B 420 (1994) 
583. 

R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and 
Schiff, Nucl. Phys. B 484 (1997) 265. 
G. Zakharov, JETP Lett. 65 (1997) 615. 
A. Wiedemann, Nucl. Phys. B 588 (2000) 303. 
Gyulassy, P. Levai and I. Vitev, Nucl. Phys. B 594 
(2001) 371. 

R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, 
Phys. Rev. C 58 (1998) 1706. 

M. Gyulassy, I. Vitev and X. N. Wang, Phys. Rev. Lett. 
86 (2001) 2537. 

X. N. Wang, Z. Huang and I. Sarcevic, Phys. Rev. Lett. 
77 (1996) 231. 
M. Gyulas sy, 



P. Levai and I. Vitev, 



th/0112071 



arXivmucl- 



R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, 
JHEP 0109 (2001) 033. 

U. A. Wiedemann, Nucl. Phys. A 690 (2001) 731. 

X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85 (2000) 

3591 

J. W. Qiu and G. Sterman, |arXiv:hep-ph/011100S . 



F. Arleo. talk at CERN TH workshop. 11-15 March 2002. 



http://wwwth.ccrn.ch/lhcworkshop/lhcworkshop01.html 



and |arXiv:hep-ph/0201066 . 
J. Binnewies, B. A. Kniehl and G. Kramer, Z. Phys. C 
65 (1995) 471. 

K. J. Eskola, K. Kajantie, P. V. Ruuskanen and 

K. Tuominen, Nucl. Phys. B 570 (2000) 379 

C.A. Salgado and U.A. Wiedemann, in preparation. 

G. David [PHENIX Collaboration], Nucl. Phys. A 698 

(2002) 227. 



